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ABSTRACT 

We study the disk and accretion properties of young stars in the NGC 2068 and NGC 
2071 clusters. Using low-resolution optical spectra, we define a membership sample and 
determine an age for the region of ~2 Myr. Using high-resolution spectra of the Ha 
line we study the accretion activity of these likely members and also examine the disk 
properties of the likely members using IRAC and MIPS mid-infrared photometry. A 
substantial fraction (79%) of the 67 members have an infrared excess while all of the 
stars with significant infrared excess show evidence for active accretion. We find three 
populations of evolved disks (IRAC-weak, MlPS-weak and transition disks) all of which 
show decreased accretion activity in addition to the evidence for evolution in the dust 
disk. 
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1. Introduction 

Circumstellar disks are a common, if not universal, part of the star formation process. The 
dust and gas within them provide the building blocks for planetary systems , hence their study is 
of critical importance for understanding planet formation processes and timescales. A considerable 
amount of work has been done in the past 20 years to inventory circumstellar disks in star forming 
regions and young stellar clusters, constraining their properties and lifetimes as a function of stellar 
age, mass, and environment. The advent of the Spitzer Space Telescope, with its unprecedented 
mid-infrared sensitivity, has opened up the field in this regard and is now providing the most 
complete census of disks to date. 

In young clusters the disk fraction is a strong function of age, decreasing fr om ~ 80% for clusters 



at 1 Myr down to almost no stars with disks for clusters older than 10 Myr (jHaisch. Lada &: Lada 
200ll ; iHernandez et al.1 120071 ). However, large variations in the disk properties exist at any age. 
Observations of the infrared spectral energy distribution (SED), which traces dust emission from 



the disk within a few AU of the star, show a wide range of excess for single clu s ters (IMegeath et al 



2005 



Sicilia-Aguilar et al 



Winston et al.ll2007 



2006 



2006 



Dahm Sz Hillenbrand 



Balo get al 



Hernandez et al 



20071 : ICurrie et al.ll2007l : iGorlova et al 



20071) as well a s a decrease in the median exces s with age (ISicilia-Aguilar et al 



20071 : lLada et al.ll2006l ). Submillimeter observations (IBeckwith Sargent 



2007 
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ture (Bouwman et al 



199ll; Mivake fe Nakagawal 19931; lAndrews k Williams 



2001 



Kessler-Silacci et al 



200 



6 



2007) and spectra o f the 10/um silicate fea- 



Sargent et al.l 120061 ) also show a large vari- 



ation in the dusts properties between circumstellar disks. The exact nature of this evolution, and 
why some stars lose their disks before others, has not been completely solved. Strong dependences 
on the steller mass have been found, particularly at ages >2-3 Myr, where les s massive stars tend 
to hold onto their disks longer (e.g. ILada et al.1 120061 : ICarpenter et al.l 120061 ) . The role of envi- 



ronment is less clear; while some disks are clearly destroyed more quickly in strong p hotoionizing 



envir onment, there is some evidence that the overall disk frequency changes little (jBalog et al 
20071 ). 



Much of the recent Spitzer studies concerning disk evolution have focused on somewhat older re- 

gions with ages >2 Myr (|Lada et al.lboOfil : IMegeath et al 1200.4 1 Isicilia-Aguilar et al.ll2006l : lDahm Sz Hillenbrand! 
20071 ) . In order to investigate earlier disk evolution, we present here an analysis of the stellar pop- 
ulation of a young embedded cluster associated with the NGC 2068 and 2071 reflection nebulae 
(hereafter treated as one region NGC 2068/71). As a site of ongoing star- formation with hundreds 
of young stellar objects, this region provide a unique environment to study early disk evolution. As 
part of the much larger Orion B molecular cloud, they are located at a distance of approximately 
400 pc ([Anthony- Twaroglll982l ). The large number of source s with infrared excess indicate it is still 



2001 



a ver y young cluster. Studies of dense gas (jLada et al.lll99ll ) and dust continuum ([Johnstone et al 



) reveal a large number of dense molecular clumps and cores, consistent with the youthful 



nature of this region. 

We have analyzed optical spectra along with optical, near and mid-infrared photometry to 
identify members of the NGC 2068/71 cluster and study the stellar and circumstellar properties 
of these members. Section [2] describes the data we have obtained, section [3] describes our results 
derived from these data, and section [J] includes a discussion of the implications of these results. 



2. Data 



Low-resolution optical spectra were taken with the Hectospec instrument on the MMT. Hec- 
tospec is a fiber- fed multi-obi ect spectrograph that uses 300 fibers to cover a 1° field of view 
(jFabricant et al.ll2005l ). These spectra have a resolution of R=1200 and cover the wavelength range 
3700-9150A. They were obtained on Nov 5,2004 and Oct 23,2005 with an exposure time of 3x900 
sec. Our sample selection methodology will be discussed below. 

High-resolution optical spectra were taken with the Hectochelle instrument on the MMT. 
Hectochelle uses 240 of Hectospec' s 300 fibers to cover a similar field of view at much higher 
resolution ([Szentgyorgyi et al.lll998l ). Spectra were centered on the Ha line at 6563A covering a 
wavelength range of 6460-6648A with a resolution of R=32,000. Two fields, each with an exposure 
time of 3x1200 sec, were observed on Dec 3,2004 and Nov 14,2005. Five stars (#177, 326, 458, 618 
and 1116) were included in both high-resolution observing runs in order to examine the variability 
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in the Hq profile of these objects. 



Near-infrared JHK S photometry was taken from the Two- Micron All Sky Survey (2MASS). 
Mid-infrared photometry was obtained using the IRAC and MIPS cameras aboard the Spitzer 
Space telescope. The IRAC observations were taken as part of GTO program 29, IRAC Survey of 
the Orion Molecular Clouds. IRAC images were taken in High-Dynamic Range (HDR) mode with 
exposures of 10.4 and 0.4 seconds ([Fazio et al.ll2004l ). Individual BCD frames (pipeline version 
S9.5.0) were mosaicked using a custom IDL program creating separate mosaics from the short 
exposu res and the long exposu res. Source selection on these mosaics was completed using Phot Vis 
vl.08 (jGutermuth et al.l 120041 ) . Photometry was derived from individual BCD frames using the 
highest signal to noise data that did not have any saturated pixels (Megeath et al. in prep). The 
photometry used an aperture of 2.4", and a sky annulus with inner and outer radii of 2.4" and 
7.2". Magnitude zero points of 21.93, 21.26, 19.08, 19.44, including the aperture correction, were 
used for these data. 

The MIPS observations were taken in scan mode using medium scan rate and half-array offsets, 
with a total areal coverage of about 40'x45' common to the three arrays. The total effective exposure 
times are 80, 40, and 8 seconds at 24, 70, and 160 microns. Since only a few Class II sources were 
detected at 70 microns and none at 160, we have omitted those data in this study. The images 
were reduced using the MIPS instrument team Data Analysis Tool, which includes basic calibration 
as well as mosaicking of individual exposures (IGordon et al J 12003 ). The final mosaic was created 



with half-pixel subsampling. We then performed photometry on individual point sources in the 
24 micron mosaic using PSF fitting with an emipirical PSF constructed from bright isolated point 
sources in the field. A PSF fit radius of 5.6" and sky annulus of 15- 22.5" were used. An ape rture 
correction of 1.73, as estimated from the STiny Tim PSF model (e.g. lEngelbracht et al.ll2007l ). was 
applied to the raw photometry. Fluxes were converted to absolute units using a conversion factor 
of 150 DN/s/mJy. 



Optical photometry was taken from the Sloan Digital Sky Survey (SDSS) ([York et al.ll2000l ). 
The Orion field was observed as part of an early calib ration run through the Galactic plane in the fall 
of 1998 and the fall of 1999 OFinkbeiner et al.ll2004l ). SDSS photometry includes measurements at 
ugriz, whose effective wavelengths are 3540, 4760, 6290, 7690, 9250A. In total we have photometry 
extending from 0.35-24/iin over a field roughly 0.8° x 1.2° centered on the NGC 2068 and 2071 
nebulae. The MIPS field is slightly smaller, which only affects a few of our selected cluster members. 
Table [T] lists the photometry for the cluster members and Figure Q] shows our cluster members 
overlayed on an optical image of the cluster. 



- 4 - 



3. Results 

3.1. Membership Selection 

Targets for the spectroscopic survey were chosen based on their position on the r vs. r-i 
color magnitude plot above the expected position of the main sequence at the distance of NGC 
2068/71, as seen in Figure [2j We restricted ourselves to stars with r magnitude between 13 and 
18.5 to exclude stars that would saturate and those that are too faint to be easily observed. The 
faint limit corresponds to an unextincted 2 Myr old M6 star at 400pc. Of these stars many are 
foreground/background stars unrelated to the cluster. The cluster members were identified from 
the presence of Lithium (6707 A) absorption and Ha emission in the low-resolution spectra. Since 
Ha emission can be seen in older stars, such as dMe stars, we use strong Li absorption as the 
primary indicator of youth, while Ha and other emission lines are used as a secondary indicator. 
Lithium absorption is an accurate indicator of youth in low-mass stars, but is more ambiguous 
for G type stars and earlier; with smaller convective envelopes these stars do not burn the Li in 
their atmospheres as quickly and Li absorption can be observed in older stars. Along with our flux 
limit, this focuses our survey to low mass stars of spectral type K to mid M. We do not attempt 
to characterize the high mass or brown dwarf population of young stars within NGC 2068/71. 

Our survey is also limited by extinction. Since we define our members based on optical spectra, 
we cannot survey the deeply embedded population. Figure [3] shows a J vs. J-H color-magnitude 
diagram, similar to Figure El with a reddening vector extending from a 3 Myr old MO star at 400 pc. 
Our sample is limited to Ay < 5 although many sources further along the reddening vector have 
excess emission at 3.6-24^m consistent with circumstellar disks and hence are probably embedded 
cluster members. A J-H vs. H-K color-color diagram for the cluster members is shown in Figured! 
The cluster members without a near-infrared excess fall within Ay = 5 of the main sequence locus 
due to extinction limit imposed by the use of optical spectra. 

Table [2] lists the stars identified as likely cluster members along with their spectral types and Li 
and Ha equivalent widths. Of the 453 stars observed spectroscopically, 69 are likely members. They 
are a mix of classical and weak T-Tauri stars (CTTS,WTTS respectively) and range in spectral 
type from G6-M6. Two are continuum stars that do not show any significant features that could 
be used to derive a spectral type. They are added as members, but are not included in the rest of 
our analysis. CTTS and WTTS are sep arated based on their Ha equivalent width, using the Ha 



EW versus spectral type boundary from lWhite Basril (120031 ). Resolved Ha profiles are needed 



to best discriminant between accretors and non-accretors, as will be discussed later. 



3.2. HR Diagram 



To construct an HR diagram and estimate an age for NGC 2068/71, the luminosity and effec- 
tive temperature are derived for each star. The effective temperature is derived from the strength 
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of multiple absorption bands and lines within the low-resolution spectra, base d on the prescription 



of Si cilia-Aguilar et al.l (|2004l ) with a few modifications. In the prescription of ISicilia-Aguilar et al 



(|2004l ) the spectra are split into broad categories of early (<K1), medium (K1-M0) and late (M0- 
M6) spectral types, and separate bands/lines are used within each broad category to derive an 
effective temperature. Indices were formed for these bands based on the ratio of flux in the feature 
to flux in the continuum. We did not use the Ti02 (4975A) band for the medium spectral types 
or the TiOl (4775A) band for late spectral types because they produced T e ff that were incon- 
sistent with the other bands. The spurious temperature from the TiOl band could be caused by 
low signal to noise in the later spectral types or veiling from accretion shocks, while the spurious 
temperature from the Ti02 band could be caused by H/3 emission in the continuum interval next 
to the feature. The Ti03 (6185A) and Ti04 (7140A) bands are best used for stars with TiOl 
index > 1.2 and Ti02 index > 1.26 indicative of later spectral types, but we used them when- 
ever possible because of their tight correlation between band strength and T e ff. When multiple 
bands were used to determine the T e ff of the star the spread in T e jf was used to estimate the 
uncertainty in spectral type. S pectral types of the m embers were visually confirmed by compar- 
ing to the standard spectra of lAllen fc Stroml (|1995l ) . The spectral types of M stars were also 
double-c hecked with the Call and CaH spectral bands of I Allen fc Stroml (119951 ). as well as the VO 



the conversion of 



bands of Iwilking etaD (hood). Effective tem peratures were transformed to spectral types using 



Kenvon fc Hartmannl (|1995l ). 



Veiling can sytematically effect our spectral typing since it will reduce the band strengths, 
resulting in an earlier spectral type. We can estimate the effect of veiling for those star where 
measur ements of the veiling e xist (see section 13. 4p . Since veiling is roughly constant redward of 
5000A (jBasri fc Batalhalll990l ) our measurements of the veiling at 6500A should be equal to the 
veiling at the bands that are used to measure spectral type. For these stars the effect of veiling is 
on the order of the spectral type uncertainty. Overall we do not account for veiling when deriving 
spectral type, but we do not expect it to substantially change our results. 

The luminosity of each star was calculated from its optic al photometry, tr ansformed from SDSS 
ugriz magnitudes to Johnsons-Cousins UBVRI magnitudes (jJordi et a l. 2006) . The extinction was 



derive d using the observed R-I color along with the int rinsic R-I color from iKenvon &: H artmann 



1995) and the extinction law of lRieke &: Lebofskyl (119851 ). We chose the extinc tion law of 



19851 ) because it forms the basis for more detailed extinction law studies (jCardelli et al 



Ric 



ce Lebofskv 



19891 ) 



It also is comparable t o recent determinations of the extinction law in the 2MASS passbands 
(jlndebetouw et al.ll2005l ) allowing us to use it when deriving reddening from R-I colors and from 
JHK colors (discussed below). T he luminosity was then dete rmined from the dereddened I band 
using bolometric corrections from IKenvon &: Hartmannl (|1995l ) and a distance of 400 pc. Two ob- 
jects, #984 and 458, have no optical photometry and no optical luminosity is calculated. Table [3] 
lists the luminosities, effective temperatures and radii, derived assuming L=47ri? 2 crT e 4 // . The typ- 
ical uncertainty in effective temperature is 5%, while the typical uncertainty in luminosity, which 
includes uncertainty in T e ff and photometric uncertainty, is 10%. 
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The luminosity can also be derived from the 2MASS J band magnitude. To determine the 
amount of extinction for each star, the observed J-H color was dereddened to the intrinsic J-H 



color, taken from IKenvon fc Hartmann (1995), for the WTTS while the C TTS were dereddened 
to the CTTS locus from iMever et al.l ([19971 ). Bolometric corrections from IKenvon fc Hartmann 
(|1995l ) were used and a distance of 400 pc was assumed. There is no systematic difference between 
the extinction or luminosity derived from the optical colors or the near-infrared colors. For three 
stars where dereddening based on R-I colors resulted in negative extinction, i.e. the observed colors 
are bluer than the intrinsic colors, the extinction was instead derived from the near-infrared. Two 
of these stars have la te spectral types (M5-M6) for which a different choice of intrinsic colors (e.g. 
Luhman et al.l (|2003l )) would result in positive extinction. For these stars the extinction is close 
to zero and a change in the intrinsic colors will have a significant effect on the reddening derived. 
Object 1171 is much brighter in the optical than in the near-infrared, possibily due to variability, 
and the extinction and luminosity are both derived from the near-infrared. For objects #984 and 
458 the extinction and luminosity were derived from the near-infrared photometry because of the 
lack of optical photometry. Object 458 appears extended in the 2MASS images giving it anamolous 
colors as seen in figure [H and its extinction is assumed to be zero even though it is likely to be 
much higher. Since the extinction is most likely non-zero the luminosity we derive is a lower limit. 
For the rest of the stars we use the optically derived luminosities for our analysis. 

Figure 5] sho ws the HR diagram for NGC 2068/71 along with theoretical isochrones from 
Siess et al.l ((2000). The members are spread from above the 1 Myr isochrone to below the 10 Myr 
isochrone, a spread beyond our errors in lum inosity and temperature and typically seen in HR 
diagrams of similar regions ( Hillenbrand! 1997 ). There appears to be no significant age difference 
between the CTTS or WTTS. The median age for the cluster is 2 Myr, with the majority of the 
stars lying between the 1 and 3 Myr isochrones. Using optical spectroscopy limits our study to 
the stars that are not heavily embedded within the molecular cloud. Our measured age of 2 Myr 
may be an overestimate of the true mean age if the substantial population of embedded stars is 
systematically younger. Without any information about their luminosities or T e ff we use an age 
of 2 Myr for the rest of our analysis. 

When studying the age and masses of the stars in NGC 2068/71 we hav e chosen to use the 
theor etical isochrones of ISiess et al.1 (120001 ) . Other isochrones, such as those by (jBaraffe et al .1 11998) 
and ([D'Antona &: Mazzitellil 119971 ) would result in a different estimate for the age of the cluster 
and the masses of indivual stars. The D'Antona & Mazzitelli tracks tend to produce younger 
ages, with most of the NGC 2068/71 members lying above the 1 Myr isochrone. The Baraffe et 
al. models agree more accu rately with dynamical mass estimates of binaries, with large scatter, 
([Hillenbrand &: Whitd [2004J) although the models do not extend to our most massive members. 
The Siess et al. tracks measure sytematically low masses, similar to the D'Antona & Mazzitelli 
tracks, although the tracks do follow the locus of cluster members in the HR diagram. We do not 
account for these systematic effects but only note that the age of the cluster and the masses of 
individual stars is model dependent. 
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Contamination from the nearby older, 4-6 Myr, OBlb association might affect our estimate 
of the age of NGC 2068/71. Without proper motion it is difficult to distinguish members of NGC 
2068/71 from members of OBlb that have strayed into our field of view. We ca n estimate the 



contam ination by considering the radial distribution of Orion OBlb members from iBriceho et al 



(|2005l ). Based on the size of our NGC 2068/71 field and distance from the center of the OBlb 
association we expect roughly one star from the OBlb association to fall within our field. This will 
not substantially affect our estimate of the age of the cluster, or our study of the disk population. 

The mass of each star was derived based on its position in the HR diagram by interpolating 
along a theoretical isochrone between the nearest model masses using effective temperature. This 
was done for the isochrone above and below each member, and the two mass estimates were aver- 
aged. Isochrones at 1, 2, 3, 5, 10, and 20 Myr were used while stars above the 1 Myr isochrone were 
assumed to lie on the 1 Myr isochrone and stars below the 20 Myr isochrone were assumed to lie on 
the 20 Myr isochrone. There is little difference in the mass estimate when using the two isochrones 
because many of these stars are young enough to still lie on vertical convective tracks. For those 
higher mass stars that lie on older isochrones there is more uncertainty in the mass estimate be- 
cause this position in the HR diagram corresponds to an evolution along constant luminosity rather 
than constant T e jf. For these stars, typically with M> 1M & , we interpolated along each isochrone 
using luminosity instead of T e ff. Since most of our stars are K to mid-M this only applied to the 
handful of stars with the earliest spectral types. The typical uncertainty in the mass estimates is 
10%, based on the uncertainty in effective temperature, ignoring systematic uncertainties in the 
theoretical models. Masses are listed in Table El 



3.3. Radial and Rotational Velocity 



Radial velocities were derived for those stars with high-resolution spectra using the cross- 
corre lation program fxcor in IRAF, along with 3 late G, early K type standard stars. lHartmann et al. 
(|1986l ) describe in detail the methodology used to derive radial and rotational velocities. Our re- 
sults are given in Table [H Two objects (#584,373) do not have any spectral lines above the noise 
that can be used to measured an accurate radial velocity, and no radial velocity is recorded. The 
distribution of radial velocities for all members is shown in Figure [61 along with just the CTTS, 
which show a similar distribution. A gaussian fit to the radial velocities of the members has a 
center at 27.1 km/sec and a standard deviation of 1.3 km/sec, which is similar to the uncertainties 
in the rad ial velocities. A simila r spread is measured in the radial velocities of Taurus and ONC 



members (lHartmann et al 



1986 



Sicilia-Aguilar et al.1 120051 ) . The radial motion of the members 



is co nsistent with the motion of the surrounding molecular cloud (jLada et al 



1991 



Wilson et al 



2005). Given the youth of the members, they have not had much time to travel outside of the cloud 
in which they were born. If we assume the velocity spread represents the difference between stellar 
and cloud velocities , it would take 5 Myr to travel 5 pc, roughly the size of the molecular cloud 
dwilson et alJbood ). 
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Rotational velocities were derived from the width of the cross-correlation function determined 
by the program fxcor. The relationship between the width of the cross-correlation function and 
the rotational velocity is determined by artificially broadening spectral standards, which all have a 
rotational velocity below the resolution limit of 8 km/sec. A quadratic function was fit to the width 
versus vsini, and the average of these fits for the three spectral standards was used to convert the 
width of the correlation function for the object spectra to vsini values. The Ha line was excluded 
when measuring the rotational velocity since it is not broadened by rotation alone. Rotational 
velocities, listed in Table [U are plotted as a function of break-up velocity (vbreakup = (GM / R) l l 2 ) 
in Figure [3 Most of the members have vsini below the resolution limit of 8 km/sec and almost all 
have a rotational velocity which is a small fraction of the break-up velocity. The two stars with 
the highest fractional rotational velocity have noisy continua, which could lead to a spuriously high 
rotational ve locity. The rotation rates are consistent with the small vsini seen in other regions of 



similar ages (jSicilia-Aguilar et al. 



2005, 



2006; 



Clarke fe Bouvierll200Ch . 



3.4. Accretion Activity 

The Ha line is a sensitive measure of accretion activity. For our sample the equivalent width 
(EW) of the Ha line, measured from the low-resolution spectra and listed in Table [21 was used 
to differentiate between actively accreting stars (CTTS) and stars that are not actively accreting 
(WTTS). We hav e used the EW boun dary as a function of spectral type between CTTS and 



WTTS defined bv lWhite Basril (|2003l ). Self-absorption generally does not affect the EW. When 
the absorption is the strongest, the emission is also very strong because the accretion and mass loss 
that the two components trace are related to each other. With this in mind resolved Ha profiles 
are still the preferred accretion diagnostic. 

High-resolution hectochelle spectra, examples of which are shown in Figure EJ were available 
for 34/67 of our likely members, as well as our two continuum stars. For these stars the Ha profile 
was used to delineate active accretors from non-accretors instead of the Ha EW. A full width at 
10% maximum flux of 250 km/sec was used to separate accretors and non-accretors. Measuring the 
width at 10% maximum flux rather than half maximum helps to avoid any self-abs orption features 



i n the profile while still measuring the profile above the noise of the continuum. IWhite fc Basri 



(120031 ) use a boundary of 270 km/sec, while we choose 250 km/sec to consistently include veiled 
stars as actively accreting, see below and Figured! Four exceptions (#1262, 326, 994, 998 )out of 
the 35 members with high-resolution spectra show signs of accretion although their high-resolution 
Ha EW lies below the CTTS boundary. Object 1262 (Ha EW=8A FW 10% =340 km/sec) has 
a very large absorption dip in the center of the Ha profile, reducing the flux in the line, and the 
measured EW, while object 326 (HaEW=5.9A, FW 10%=170 km/sec) has an inverse P Cygni 
profile. The presence of a strong absorption feature, such as the inverse P Cygni profiles in objects 
326 and 515 in figure [8l can lower the velocity width below our threshold of 250 km/sec. In these 
cases the presence of strong absorption features due to accretion flows, as well as an asymmetric 
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profile, indicate that accretion is still ongoing. Object 984 has a very noisy high-resolution spectra, 
possibly due to its high reddening, Ay = 4.5, making it difficult to evaluate its accretion status. 
Object 998 also has a noisy profile and the high-resolution spectra of these two stars was binned to 
increase the signal to noise. The width of the Ha profile is above our threshold for active accretion; 
however the accretion status of these stars is still questionable without a higher signal to noise Ha 
profile, or another sign of active accretion such as a U band excess. Unfortunately the red leak in 
the SDSS u-band prevents us from accuratly measuring the U-band excess in this filter. The Ha 
EWs measured from the high-resolution spectra, listed in Table [U are consistent with the Ha EW 
measured from the low-resolution spectra. 

Continuum veiling of the photospheric absorption features provides an additional measure of 
accretion activity, albeit one that is less sensitive than the Ha profile. Accretion shocks produce 
a featureless continuum that adds to the photosphere, reducing the strength of photospheric lines 
relative to the observed continuum. The change in these photospheric features can be measured 
b y comparing a veiled s pectra with an unveiled spectra of the same spectral type. The method 
of lHartigan et al.1 (|1989l ) was used to calculate the veiling (r = F excess / -F cori i) over the wavelength 
range 6490-6501A. This wavelength range was chosen because it contains a number of strong stellar 
lines while other wavelength ranges in our high-resolution spectra did not consistently have strong 
lines useful for veiling and were ignored. We added a constant flux to the unveiled spectra until 
the strength of the photospheric features relative to the continuum matches that of the veiled 
spectra. This constant flux level relative to the continuum is designated r. The uncertainty in 
the veiling is approximately 0.15, and was estimated by comparing standards of similar spectral 
types to each other. Standards of spectral type K1.5,K2.5,K7,M0,M1 were used and veiled spectra 
were compared to standards within one spectral type. For the veiled stars that have a measurable 
rotational velocity, the standard was artificially broadened to match the rotational velocity of the 
vei led star. We found ve iling measures ranging from 0.31 to 3.2, similar to those found near 6500A 
bv lHartigan et al.l (|l99lh . 



Accretion rates were estimated from the veiling measures using the procedure of lWhite &: Hillenbrand 



(J2004J) . The conversion starts by using the R band magnitude as a measure of the continuum emis- 
sion in this wavelength range, then from the definition of r we determine the excess flux, converting 
to an accretion luminosity by assuming the total accretion flux was 11 times larger than the flux 
observed over the R band and then assuming the accretion infall comes from 3ii* to convert accre- 
tion luminosity to accretion rate. Correcting the observed accretion flux for the total accretion flux 
depends on the model of the accretion shock chosen. We use a correct ion factor of 11, a loga r ithmi c 
average of various models, which is consistent with the analysis of IWhite Hillenbrand! (|2004l ). 
The radius 3-R* is typical for the truncation of the disk by the magnetosphere, where materia l 
from the disk flows onto the magnetic field lines (jKenvon et al.lll996l ; IWhite &: Hillenbrand! 12004 ). 
Instead of using theoretical R band magnitudes, we used observed R band magnitudes, correcting 
for the contribution from the veiling and extinction. This avoids the assumption that all of the 
stars have exactly the same age, and is not subject to the systematic uncertainties associated with 
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theoretical isochrones. For those stars with no measurable veiling but signs of accretion, such as a 
wide and/or asymmetric Ha profile, an upper limit on the accretion rate was calculated assuming 
a veiling of r=0.15. Accretion rates are listed in Tableland veiling, excess flux from accretion in 
the R band, accretion luminosity and accretion rate are plotted in Figu re | 10l as a function of mass. 
The M values are typical of the range in young regions such as Taurus (jGullbring et al.lll998l ). We 
do not cover a large enough range in mass to look for any statistically significant trends. 

Five of our likely members, objects 326, 177, 458, 618 and 1116 have two high-resolution 
spectra taken a year apart. Object 618 has the largest variation with both a change in the shape 
and strength of the profile, even though the width of the profile does not significantly change, as 
seen in Figure [8l The other stars do not vary significantly in the shape or strength of the Ha 
profile. The 10% F W vary by less than 4 km/sec for all of these s tars, similar to other stars with 
similar Ha profiles (jJohns &; Basrilll995l : iJayawardhana et al.ll2006l ). 



Accretion is often associated with outflows and jets driving material away from the star 
(Kon igl fc Pudrita 120001 ) . The NGC 2068/71 cluster has a number of strong molecular outflows 



((White fc Phillips 



Haro objects (jZhao et al 



1981 



1999 



Snell fc E dwards! 1982: iBallvl Il982: iFukui et alJll98fl ) as well as Herbig 



Reipurthl Il999i ) throughout the cloud. None of these outflows 



appear to be associated with the pre-main sequence stars studied here, and may instead be related 
to the younger, less evolved stellar population of NGC 2068/71. In the optical spectra, a number 
of stars exhibit [OI] 6300A emission associated with outflows. These stars, marked in table all 
have an infrared excess and are actively accreting supporting the view that accretion and outflows 
are physically connected. 



3.5. Disk emission 



The shape of the infrared SED contains information about warm dust surroundin g the young 
star and can be used to distin guish between stars with disks and stars without disks. lLada et al 



(|2006l ) and iHernandez et al.l (j2007l ) used the slope of the SED in the IRAC bands, ajRAC, to 
distinguish between stars with strong infrared excess coming from optically thick accretion disks, 
stars with a weaker excess indicative of evolved, or "anemic", disks, and stars with no measurable 
excess associated with a lack of disk material within 1 AU. We refer to the evolved objects as 
IRAC-weak because they may have strong excess at longer wavelengths indicative of optically 
thick outer disks . Figure Qj] shows oirac versus effective temperature for NGC 2068/71. As in 
Lada et al.l (|2006l ) we find two distinct populations, one around Q-irac = —1-0 and a second around 
0(irac = —2.8. We have not dereddened t he SEDs since the low extinction (A v < 5) and relatively 
flat extinction law across the IRAC bands ([Flaherty et al.ll2007l ) means that dereddening will have a 
small effect on the SED slopes. Previously a boundary of a[RAC=-\ - 8 has been used to distinguish 
evolved disks from strong disks (|Lada et alj 120061 : IHernandez et al.l 120071 ). We use a boundary of 
o/i?AC = -l-6 because of the stars near ocirac=-1.8 that likely have some evolution in the structure 
of their disks. For our analysis stars with amAC> —1-6 were treated as having strong disks, stars 
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with —2.5 <aiRAC< — 1-6 have IRAC-weak disks and stars without disks have ccirac< —2.5. 

An IRAC color-color diagram, Figure [T2l shows a similar progression in disk properties. There 
is a population of sources with no infrared excess, a population with strong infrared excess, and 
population in between with weak but significant infrared excess. Including the 24/im band of MIPS 
adds additional information on the dust properties of the disk. Figure [131 shows the [8.0]-[24] vs. 
[4.5]-[5.8] color-color diagram. All of the detections at 24/im, 33/34 strong disks and 8/12 IRAC- 
weak disks, have some sort of disk emission as expected given the 24/tm sensitivity. Of the stars 
with an infrared excess not detected at 24/im one is outside the MIPS field of view and the other 
three are in regions of strong nebulous emission where the sensitivity is decreased. 

There is a population of sources with [8.0]-[24]<2.4, which we refer to as MlPS-weak. These 
sources are singled out because this [8.0]-[24] color indicates spectral slopes < —1, similar to \F\ oc 
A -4 / 3 , and because of the evidence for lower accretion activity, discussed below. A perfectly flat 
disk has an SED with the form \F\ oc A~ 4 / 3 , which represents the limit of a disk with a substantial 
amount of grain settling. Figure [Til shows olirac y s- chmips > the SED slopes from 3 — 8/im and 
8 — 24/im respectively. The MlPS-weak sources all have a strong excess at the IRAC bands, even 
though they have a weak excess at 24/im. The lack of MIPS weak sources with small olirac may 
be caused by the limited sensitivity at 24/im. The majority of sources show strong excesses at both 
the IRAC bands and at 24/im, while the IRAC-weak sources have a range of excesses at 24/im. 
Figure [15] shows sample SEDs for each of the types of circumstellar disks. The biggest change in 
the SED between the MIPS weak star and the strong disk star is between 8 and 24/im with the 
MlPS-weak disk having a much steeper SED beyond 6/im. The SED of the IRAC-weak sources 
looks similar to the SED of the strong disk but with a steeper 2-8/im slope. 

Stars contaminated by background PAH emission may appear to have an excess at [8.0], even 
though they have no circumstellar material within a few AU. Three stars, #802, 984, 1078, are in 
regions of high nebular emission and have an excess primarily at [8.0]. Excess emission at 24/im 
would help to better characterize the presence of a disk, but none of these stars are detected at [24] 
possibly due to the high background emission surrounding the stars. We leave them as IRAC-weak 
disks but note that their disk status is still uncertain 

Including all stars with olirac> —2.5 as having a disk as well as the stars with cxirac< —2.5 
but a 24/im excess indicative of an optically thick outer disk, 79% (53/67) of the spectroscopically 
confirmed members have disks. This disk fraction may not accurately reflect the disk fraction in 
the entire cluster, since sources with interesting mid-infrared colors were added to our spectroscopic 
sample. To constrain the true disk fraction we consider the photometry of those confirmed members, 
as well as stars with similar colors. In the r vs. r-i color-magnitude diagram (shown in Figure [2]), 
there are 397 stars that sit above the main sequence locus and have r< 18.6, 230 of which have 
been observed spectroscopically. Of these 230, 44 are members with disks (ajRAc> —2.5) and 14 
are members without disks. Our sample of members without disks is unbiased since we have no 
prior information on whether or not they are members. Scaling up from the spectroscopic sample 
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to the photometric sample, we expect 24 stars without disks in the entire sample of 397 stars. If we 
assume that we have found all the members with disks, then out of the entire sample there would 
be 24 stars without a disk, and 44 stars with a disk for a disk fraction of 44/68 (65%). We have 
likely not found all the members with disks, and the true disk fraction is likely higher than 65%. 
Out of the 397 stars, 58 have infrared colors consistent with disk emission. If we instead assume 
that we have found all of the members without disks and that all of the objects with an infrared 
excess are members, then the disk fraction would be 58/72 (81%). Some of the objects with an 
infrared excess may be background AGN or AGB stars and it is also unlikely that we have found 
all of the members without disks so we expect the true disk fraction to be less than 81%. This does 
not count stars that could be missing material within 1 AU but the small number of these objects 
(see § I3.5.ip will not have a significant effect. The spectroscopic sample has a disk fraction of 79%, 
while we expect the true fraction for the entire cluster to be in the range 65-81%. 

Another e stima te of the disk fraction can be made from the X-ray emitting stars observed by 



Skinner et al.1 (120071 ). X-ray emission as a tracer of young stars is not biased towards stars with 
disks. We detect 23 of the 33 X-ray emitting stars with IRAC, mainly at [3.6], [4.5]. We use these 
two colors, with a boundary of [3.6]-[4.5]=0.2 to differentiate stars with or without disks. One star 
is not detected by 2MASS, which combined with its faint IRAC flux suggests it is an extragalactic 
source and we exclude it. Where possible we use detections in other bands to better characterize 
the presence of a disk around these stars. Two stars with [3.6]- [4.5] > 0.2 have colors consistent 
with highly reddened photospheres while two stars with [3.6]- [4.5] < 0.2, which are both part of our 
spectroscpoic sample (#984,1078), have a slight excess at [8.0]. In total there are 14/22 (64 ± 17%) 
stars with an infrared excess in this field. This is consistent with our previous estimate of the true 
disk fraction. Ten of these stars, six of which have an infrared excess, are deeply embedded and 
are not detected in the optical. The X-ray fields are in areas of high background emission at 24/iin 
reducing the sensitivity at this band and making it difficult to estimate the fraction of evolved disks 
in this field. 

If the deeply embedded population has a higher disk fraction than our spectroscopic sample 
then this will introduce an additional bias in our measure of the cluster disk fraction. Ices are 



expected to form on dust grains within the molecular cloud above A v =3 (IWhittet et al.ll200ll ). We 
might expect that the stars with A v >3 differ from the less embedded population because of the 
change in the surrounding molecular cloud. Of the stars with A v > 3, 93% (28/30) have disks 
while 71% (20/28) of the stars with A v < 3 have disks. This only includes stars earlier than M4; 
embedded stars later than M4 are too faint in the optical to be observed here. There does not 
appear to be a significant age difference between these two groups although the change in the disk 
fraction is suggestive of a difference between the embedded population and the revealed population. 
Observations of more heavily extincted cluster members are needed to characterize this difference 
fully. 
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3.5.1. Disk Hole Sources 



A number of stars have been observe d towards other young clusters that appear to be missin g 



dust within a few AU of the central star (| Guillen et al 



2004 



Calvet et al. 



2002 



Rice et al J 12003 ). 



These objects, often referred to as transition disks, lack infrared emission shortward of lO^m but 
exhibit strong emission longwards of 10/im. NGC 2068/71 has two such objects, objects 177 and 
281. They have [8.0]-[24]> 4 and [3.6]-[4.5]<0.1. Our high-resolution spectra allows a search for any 
weak accretion. Object 177 (FW 10%=124 km/sec, Ha EW=1.2A) does not appear to be actively 
accreting, while object 281 has a broad, asymmetric Ha profile indicating active accretion (FW 
10%=400 km/sec, Ha EW=18.lA). A veiling of r=0.36 was measured for object 281 corresponding 
to an accretion rate of M = 2.4 x 10 -9 M^yr" 1 . This accretion rate is the smallest rate measured 
among the stars in this cluster, although our measured accretion rates are biased toward those stars 
whose accretion rate is high enough to have a measurable veiling. 



3.5.2. Infrared Excess and Accretion 

The stars with strong disks (32/67, 48%), with spectral types ranging from G6-M5, all have 
signs of active accretion. They have strong Ha EW (5-300A) and broad asymmetric Ha profiles. 
The full width at 10% maximum flux for the stars with strong disks ranges from 260 to 640 
km/sec. For those stars with measurable veiling the accretion rates are in the range 7.9 x 10 -9 — 
8.6 x W~ s M^yr -1 . The stars with no disks (14/67, 21%) show no sign of active accretion, either 
in Ha EW (1.3-9.5A) or the Ha profile (FW at 10%=90-200 km/sec). Figure H shows the Ha EW 
versus FW at 10% where the strong disks and stars with no disks form two distinct groups. 

The MlPS-weak sources (8/67, 12%), with spectral types ranging from K4-M4, show signs 
of lowered accretion activity in addition to lower infrared excess. In the color-color diagrams, 
Figures [T31 [T4"l the size of the symbols scales with the Ha EW and the MlPS-weak sources have 
smaller Ha EW than the strong disks. These eight stars have Ha EW between 3 and 48A. Three 
have high-resolution spectra with 10% full-widths of 480, 340 and 450 km/sec, and accretion upper 
limits of 2.4 x 10~ 8 M Q yr~ 1 , 1.79 x lO -8 M yr -1 and 1.8 x l$- 7 M Q yr~ l . In Figure [9] some of the 
MlPS-weak sources with high-resolution spectra are distinct from the strong disks and the stars 
with no disk. The HaEW is lower than the strong disks with comparable FW at 10%. 

The IRAC-weak sources (12/67, 18%, not including #281), with sp ectral types rangin g from 



K1-M4, have -2.5 <a IRA c< -1-6 and also tend to have [8.0]-[24]> 2. lLada et al.l (|2006ft find a 
similar spread in [8.0]- [24] color for their optically thin/" anemic" disks. Of these twelve stars, six 
have high-resolution spectra that show 10% full widths for Ha of 280, 410, 280, 260, 170 and 520 
km/sec, close to the boundary between accretors and non-accretors. For these six objects only 
upper limits of 7.11,0.94,1.78,5.4,0.47,1.06 x lO _8 M yr _1 can be placed on the accretion rate 
since no veiling is measured. Figure [9] indicates they have the same accretion activity as the three 
MlPS-weak stars, separate from the strong disks and the stars with no disk. These stars have Ha 
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EW between 2 and 32A, with one star having an Ha EW of 75A. With the exception of the one 
active accretor, these stars show less accretion activity than the strong disks. 



4. Discussion 



Many studies now show that the disk fraction d ecreases rapidly with age, wi t h nearly all stars 



losing their primordial disks after roughly 10 Myr (jHaisch. Lada Lada 



2001 



Hernandez et al 



20071 ). For NGC 2068/71 we measure an age of 2 ± 1.5 Myr and a disk fraction of 79%. Our disk 
fraction is consistent with ground-ba sed estimates from JHKL fo r regions such as the Trapezium 



and NGC 2024 with ages of 1-2 Myr (jHaisch. Lada fe Lada 



2001 



In NGC 2068/71 the members with disks do not show a systematically different rotational 
velocity than those stars without disks. Rotational velocities, vsini, relative to the break-up speed 
are shown in Figure with the stars with and without disks marked as red and blue respectively. 
There appears to be no significant separation, although the substantial number of upper limits 
makes it difficult to evaluate any difference fully. It has been suggested that the slow rotation rates 
for pre-main sequence stars are due to disk-locking, where the star is connected to t he disk by its 



magn etic field lines and angular momentum is transferred from the star to the disk. iRebull et al 



(|2004l ) suggest that even if the slow rotation speed was due to disk locking, a significant difference 
between the rotation velocity of stars with disks and stars without disks would be difficult to 
measure. More data are needed to test these hypotheses in detail with the members of this region. 

The high-resolution spectra not only allow us to measure rotational velocities, but also obtain 
a more sensitive measure of the accretion activity with the Ha profile. All of the WTTS, classified 
based on Ha EW, that have a disk and high resolution spectra exhibit Ha profiles consistent 
with active accretion, except the transition disk #177. Objects #849, 458, 581 and 998 all have 
low-resolution Ha EW which would make them appear as WTTS, while their high-resolution Ha 
profiles indicate that they are still accreting material. This is especially important when examining 
the accretion activity of the evolved disks which tend to have low Ha EW although they are still 
accreting (see Figure [9]) . 

The location of our s elected members does not s trictly overlap with e ither the dense gas 
emission (jLada et al.lll99ll ) or the dense dust emission (j Johnstone et al.ll200ll ). Star formation is 



ongoing in these regions of dense gas and dust, based on observable outflows and class I sources. To 
compare the distribution of members to the dust cloud we created an extinction map (Figure [To]) 
by selecting background stars, assuming they have the same intrinsic H-K s color, and calculating 
E{H — K s ) = H — K. q — (H — K. q )n. A 30" gaussian kernel was used to smooth the map and 
the extinction law of iRieke &: Lebofskyl ()1985l ) was used to c o nvert from E(H-K S ) to Ay. A more 



detailed description of this method is given in iGoetz et al.l (120031 ). The extinction map closely 
matches the CS e mission measur e d by Lada et al.l ()199ll ). It also agrees with the lower resolution 



extinction map of iDobashi et al.1 (|2005l ). although their map does show a more extended peak at 
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"2000 = 05 47 m 36 s , <52ooo = +00 o 18'00" which could be due to our small number of background 
stars in this area. As seen in Figure [T6l our selected members follow the general shape of the 
diffuse dust trace d by the extinction ma p, but are not confined to the peaks of the dust, where the 
SCUBA clumps (I Johnstone et al.ll200ll ) and class I sources (Muzerolle et al. in prep) are located. 
Our sample may be biased toward those areas of lower extinction away from the dense gas and 
dust because the use of optical spectra for spectral type classification limits our survey to those 
young stars that are not deeply embedded within the cloud. Another possibility is that the star 
formation in this region has progressed from west to east. Stellar winds and photoionizing flux 
from the OB associations to the west of NGC 2068/71 may have been interacting with the cloud, 
causing stars to form on the western side of the cloud first and on the eastern side of the cloud 
later. A third possibility is that the selected members may have had time to disperse from their 
original birthplace and have spread throughout the cloud. For a typical velocity of 1-2 km/sec 
relative to the cloud, these members could have traveled 2-4 pc in their lifetime of 2 Myr, roughly 
the east-west extent of the members. 



Binarity may have a significant impact on disk properties, however, our d ata lack the spatial 
resolu tion to address this. One star (object #739) was observed to be a binary bv lPadgett. Strom Ghez 
(jl997l ). This binary has a separation of 0.97", which corresponds to 390 AU at a distance of 400pc. 
Interestingl y this star has a MlPS-weak di sk, which might be explained by tidal truncation of the 
outer disk. iPadgett. Strom k, Ghezl (jl997l ) estimate a binary fraction of 15%, which would imply 
that 10 of our sources are binaries. None of our stars with high-resolution spectra appear to be 
spectroscopic binaries but more epochs are needed to rule out single-line spectroscopic binaries. 



4.1. Evolved Disks 



There are a number of possible explanations for the SEDs of strong disks vs. MIPS-weak/IRAC- 
weak disks. Grain growth and settling is one process that may have a measurable effect on the 
SED of an accretion disk. Models suggest that thi s can be a very fast process p roducing noticeable 
differences in the SED in less than a million years (jDullemond &: Dominikll2005l ). The overall effect 
is to quickly deplete grains out of the upper layers of the disk, while grains in the inner disk will 
grow and settle more quickly than in the outer disk. 

As the degree of settling increases, the most direct effect is a decrease in the irradiation 
surface of the disk, which leads to less emerging flux. The SED approaches that of a perfectly 
flat disk, \F\ oc A~ 4 / 3 at wavelengths corresponding to the disk radii that are affected by settling 
(jD'Alessio et al.ll2006l ). Assuming the amount of grain growth and settling is constant throughout 
the disk, the change in flux is most noticeable for A > 10/im because these wavelengths trace regions 
in the disk where flaring becomes prominent and the surface height is most sensitive to grain growth 
and settling. The emission at A < 10/um is dominated by the wall at the dust destruction radius, 
which is less affected by the degree of flaring of the disk. For a substantial amount of grain growth 
and settling the height of the wall will decrease by a factor of 2, while the height of the outer 
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disk can decrease by more than an order of magnitude (jD'Alessio et al.l 120061 ) . Theory predicts 
grain growth and settling should occur more rapidly in the in ner disk than the outer disk since the 
orbital period is shorter (jDullemond Dominiki 12004 . 12005). The IRAC-weak stars may have a 
substantial amount of grain growth and settling in the inner disk, whereas the MlPS-weak sources 
may have advanced settling over a wider range of radii. 

The predicted SEDs for grain growth and settling are not unlike those corresponding to a 
decrease in the accretion rate. A lowered accretion rate will affect the disk in three ways: (1) 
decrease the surface density, (2) decrease the viscous heating and (3) decrease the accretion lumi- 
nosity emitted by the stellar accretion shock. A decrease in surface density will lower the flux at 
all wavelengths, while the smaller viscous heating will lead to a cooler midplane temperature and 
a change in the scale height of the disk similar to that of grain growth and settling. The decrease 
in the accretion luminosity will have an effect on the inner disk edge which dominates the flux at 
A < Wfim. Accretion luminosity, along with stellar irradiation, illuminates the inner dust edge of 
the disk. A decrease in accretion luminosity will cause the dust destruction radius to move closer to 
the star. While the inner edge will still be at 1500K, its solid angle will decrease and the emission 
from the inner edge will decrease. For M > 10~ 8 M^yr^ 1 the accretion luminosity is comparable 
to the stellar luminosity and a change in accretion rate leads to a change in the radius of the wall, 
and a change in the flux in the IRAC bands. Since the accretion luminosity {L = GMM / R) and 
the stellar luminosity roughly scale with stellar mass, the dependence on accretion rate dominates 
the ratio of accretion to stellar luminosity. For M < 10 -8 M & yr~ l the change in the radius of the 
wall with accretion rate is smaller because the stellar luminosity dominates. 

A change in the SED slope could be due to grain growth and settling or a decreased accretion 
rate separately or it could be caused by a combination of the two. As the amount of settling inc reases 
the effect of a change in the accretion rate becomes more significant (jD'Alessio et al.ll2006l ). For 
high accretion rates, viscous heating dominates over other forms of heating in the surface layers. 
Since viscous heating has little dependence on the grains in these layers, there is little change in 
the SED when grains start to grow. For lower accretion rates, stellar irradiation dominates the 
heating of the surface layers, and this form of heating is very sensitive to the grain distribution. The 
MlPS-weak and IRAC-weak sources may be separating themselves from the strong disks because 
their low accretion rates lead to a more significant difference in the SED. However, we need more 
estimates of M to directly compare to disk models. 

This picture is confused by the 10/xm silicate feature which overlaps with the 8//m band of 
IRAC. A strong silicate emission feature will artificially increase ojrac despite a steeper continuum 
slope. As the accretion rate decreases the intensity of this feature increases due to the contrast in 
temperature between the surface (heated by stell ar irradiation) and the midplane (heated by viscous 
dissipation and emission from the upper layers) (jD'Alessio et al.l l2006). The silicate feature is most 
sensitive to grain growth in the upper layers of the disk and does not trace grain growth in the 
midplane. Spectroscopy of the Wfim silicate feature wo uld help to break the degeneracy betwee n 
grain growth and settling and a change in accretion rate (jD'Alessio et al.ll2006l ; iFurlan et al.ll2006l ). 
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as well as seperate out its contibution to (Mirac- 



Another possibility for the observed SEDs of the MlPS-weak stars is the truncation of the 
disk by a close binary. A close binary will truncate the disk around the primary at 1/2 to 1/3 
the semimajor axis of the binary system (jArtymowicz Lub owl 1 19941 ). The majority of the 24/im 
emiss ion comes from within 0.5-5 AU, depending on the spectral type and degree of settling in the 
disk (jD'Alessio et all 120061 ) . A binary companion separated by ~1-15AU would truncate the disk 
and reduce the 24/im emission without affecting the shorter wavelength excess. The frequency of 
MlPS-wea k stars in NGC 2068/71 is roughly consistent with the frequency of binaries at separations 
~1-15AU (jMathieul Il994l ). The smaller resevoir of material available to accrete onto the central 
star may lead the accretion rate to decrease earlier than it would for stars with disks that have 
not been truncated. The smaller disks imply a smaller viscous timescale for the evolution of the 
disk, which may lead it to be quickly dissipated unless it is ressuplied by a circumbinary disk. The 
binary fraction within a cluster seems to be a function of environment, rather than of age, with 
low-desnity clusters such as T aurus and Ophiuchu s showing a higher binary fraction than dense 
clusters such as the Trapezium (jMathieu et al.ll2000l ). Millimeter observations find that disk masses 
are much lower for bi naries separated by 1-100 AU than single stars or wide binaries, suggesting 
truncation of the disk (jMathieu et al.ll2000l ). If MlPS-weak disks are due to the presence of a binary 
then we may expect the MlPS-weak fraction to be a strong function of environment rather than 
age. Interestingly one of our MlPS-weak stars (#739) is a known binary ([Padgett. Strom fc Ghez 
19971 ). although its separation (390AU) is likely to be too large for any truncation of the disk to 
affect the 24/im flux. 

The third group of evolved disks are the transition disks, with AU-scale optically thin or 
evacuated holes in the center of the disks. We have two such objects, #177, and 281, making 
up 3% of our spectroscopic s ample . One appears to be actively accreting (#281) while the other 
is not. ISicilia-Aguilar et al.l (|2006l ). who define transition disks as having no excess shortward of 



6/mi, find that half of the transition disks in Tr 37 are still actively accreting, based on high- 
resolution H a profiles. The transition disks around TW Hya, GM Aur and DM Ta u are all seen to 



be ac creting (jCalvet et al.ll2002l . 120051 ) while CoKu Tau/4 is not actively accreting (jD'Alessio et al 
20051 ). A ccretion rates of 10~ 8 — M^yr -1 have been measured fo r GM Aur, DM Tau, and 



TW Hya (IGullbring et al.lll998l : ICalvet et alJl2005l ; lMuzerolle et alJl200d ). similar to that measured 
for object #281. 

A hole could be swept out by a planet, be caused by substantial grain growth in the inner disk 
or be due to photoevaporation of the disk. The formation of a planet would cause material outside 
the planet's orbit to gain orbital angular momentum and be pushed outwards, whi le material withi n 
the planets orbit would lose angular momentum and quickly accrete onto the star (jRice et al.ll2003l ). 
The mass of the planet and the disk will affect whether a gap ca n be created, and how the planet 
reacts to the material in the outer disk (IPapaloizou et al.1 12 007) . A planet has been invoked to 
explain the ob served inner holes in GM Aur jRice et allbooal ). CoKu Tau/4 dQuillen et alJbooj ) 
and TW Hya (jCalvet et al.l 12002) . For object # 281, the material within the planet's radius may 
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not have completely accret ed onto the star, or the p lanet may be small enough that some gas may 
be able to stream past it (|Lubow k, D'Angeldl2006l ). This gas would transport small grains into 
the gap, creating an optically thin disk, which would not be detectable from photometry alone 
(jRice et al.l 12006). It is also possible that grains within 10 AU have grown to millimeter sizes 
or larger, while the grains at larger radii have not grown as much. The growth of grains should 
proceed on an orbital timescale, which would r esult in more grain growth in the inner disk than the 



outer disk (Dullemond &z Dominik 



2004 



20051 ). although this process should not produce a sharp 



boundary in the dust disk. Any accretion onto the central star may also transport small grains 
back into this gap. A third pos sibility arises from a photoevaporative wi nd from the disk driven by 



stellar UV and X-ray photons (jClarke et al 



2001 



Alexander et al.ll20 06). When the accretion rate 



is low the wind will counteract the inward flow of material outside of the gravitational radius and 
will prevent additional gas from being transported into the inner disk. With the inner disk starved 
for material any gas will quickly accrete onto the star and disappear producing an inner hole for 
a short time before the entire disk dissipates. This process requires very low (< 10~ 10 M^yr^ 1 ) 
accretion rates in order for the wind to counteract the inward flow of material due to accretion. 
This would rule out the photoevaporative model for object #281, but it is still possible for object 
#177. 

None of these three groups of evolved disks appear to be systematically older than the rest of 
the cluster, according to the HR diagram shown in Figure Although the two transition disks 
are below the 3 Myr isochrone their small number make any statistically significant conclusion 
impossible. While the number of sources in each group is small enough that we cannot do a 
detailed study of the ages, we can say that it is unlikely that either the IRAC-weak, MlPS-weak 
or transition disks are systematically older or younger than the other sta rs in the cluster. Th is 



indicates that the disks around different stars evolve at different rates. Dullemond fc Dominik 



(|2005l ) find that grain growth proceeds much more quickly than the typical age of a T-Tauri star 
and small grains must be replenished in the disk in order to retain a significant IR excess. If 
replenishment occurs then the typical grain size in a disk does not reflect the age of the system. 
We would then not expect these stars to appear systematically older if the change in SED for the 
IRAC-weak and MlPS-weak disks is due in part to grain growth. 



IRAC-weak disks ar e observed in cluster regio ns from Taurus at 1 Myr (ILada et al.1 12006) 
to NGC 2362 at 5 Myr (|Dahm fc Hillenbrand! l2007h. The fraction of IRAC-weak disks depends 
on the spectral type range under consideration. ILada et al.1 (|2006l ) find the fraction of stars with 
IRAC-weak disks ranges from 8-25% depending on the spectral type. Later spectral types have 
a larger fraction of IRAC-weak disks although this trend may be affected by larger photometric 
uncertainties at later s pectral types blurring the boundary between IRAC-weak and diskless stars 
([Hernandez et al.ll2007l ). We examine the fraction of disked stars that have an IRAC-weak SED for 
broad spectral type ranges to look for trends with age. Here we use o.irac = -1-8 as the boundary 
for IRA C-weak stars to be co nsistent with previous studies. We consider K0-M1 stars in Ta urus at 
1 Myr (|Hartmann et alJboOfil ). NGC 2068/71 at 2 Myr, IC 348 at 2-3 Myr (ILada et alJhoOfth . a Ori 
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at 3 Myr (Hernandez et al 



2007), Tr 37 at 4 Myr (jSicilia-Aguilar et aUEOOfl ) and NGC 2362 at 5 



Myr (jDahm &: Hillenbrand! 120071 ). The fractions of disked stars with IRAC-weak SEDs are listed 
in Table [5j The fractions of cluster members with disks in these clusters decrease from 80% to 
~ 20% with age, but the fractions of disks with IRAC-weak SEDs do not show an obvious decrease 
with age, although the small number statistics for some of these clusters make this trend uncertain 



Typical Ha EW for t he IRAC-weak disks in Tr 37 and IC 348 are < 15 A (ISicilia-Aguilar et al 



2006 



Lada et al.ll2006l ). consistent with the IRAC-weak stars in NGC 2068/71, and ava i lable high- 



resolution Ha profiles indicate they are actively accreting. However iHernandez et al.l (120071 ) find 
that the median SED slope for all disks in a given region decreases with age from 1-3 Myr, which 
may contribute to the fraction of disks with IRAC-weak SEDs if the entire distribution of infrared 
excess is shifting to a smaller excess. 



Of the 54 stars with disks detected by MIPS in Tr 37, 7 have MlPS-weak disks (ISicilia-Aguilar et al 



20061 ). while out of t he 38 K0-M4 stars with disks in IC 348, 3 have MlPS-weak disks, none of which 
are earlier than Ml (jLada et alJl2006h. In a Ori. three out of 16 K0-M1 stars with disks detected 
by MIPS have MlPS-weak disks (IHernandez et al.ll2007l ). The sensitivity limit at 24/xm makes any 
trends with age harder to evaluate although we note that the fraction of disks with MlPS-weak 
SEDS does not significantly decrease with age. The Ha EW of the MlPS-weak stars in Tr 37 are 



between 4 and 47A, with most below 15A ( Sicilia-Aguilar et al. 2006). O f the stars with measured 
accretion rates, most are below lO _8 M0?/r _1 ( Sicilia-Aguilar et al.l 1200 6). This is consistent with 
the properties of MlPS-weak disks in NGC 2068/71, although overall accretion rates are lower in 
Tr 37 and the low accretion rates for the MlPS-weak stars may simply reflect viscous evolution of 
the disks. 

If disk evolution is not solely dependent on age then it may depend on environment or the 
initial conditions of the protostar. Material infalling onto the disk from the surrounding cloud could 
replenish the small grains in the disk as well as maintain the presence of a disk. The longevity of a 
disk would then depend on the resevoir of material surrounding the star. During the early stages of 
star formation material accretes from an envelope onto the disk, and this infall rate depends on the 
properties of the envelope. More material accreting onto the disk will cr eate a more massive disk, 
which can affect its subsequent evolution (IDullemond fe Dominikl 120051 ) . Unfortunately sub-mm 
observations are not sensitiv e to circumstellar disks in NGC 2068/71and cannot be used to estimate 
disk masses for these stars (1 Johnstone et al.ll200ll ). The initial an gular momentum of the system 
may influence the initial size of the disk and how quickly it evolves (jDullemond et al.ll2006l ) as well 
as affect the possibility of forming a binary, which could enhance the subsequent evolution of the 
disk. A nearby O/B star would increase the ionizing flux impinging on the surface of the disk, and 
possibly increase the speed with which the disk dissipates, although the evolved disks we observe 
do not appear preferentially close to the B stars in this cluster and with the earliest member being 
a B2 star the UV flux is not especially strong. 
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5. Conclusion 

We have studied the young stellar population of NGC 2068/71 and identified 69 members, for 
67 of which we have derived spectral type and luminosity. Our cluster members range from G6-M6 
with a median age of 2 Myr. A large fraction (79%) of these stars have an infrared excess, although 
this is likely an overestimate of the true disk fraction. Of the stars with an infrared excess and a 
high-resolution Ha profile, sensitive to low-accretion rates, all but one are actively accreting. The 
stellar members are not confined to regions of dense gas and dust, although there is some evidence 
that the disk fraction is higher in dense {Ay > 3) regions. 

Combining accretion information from the Ha line with the 3-24/xm infrared SED we find three 
populations of evolved disks: IRAC-weak, MlPS-weak and transition disks. All of the populations 
show a change in the shape of the infrared SED as well as a decrease in the accretion activity. For 
the IRAC-weak disks this could solely be due to a decrease in the accretion rate (which affects the 
structure of the disk), substantial grain growth and settling, or a combination of the two. These 
two processes could also explain the MlPS-weak SEDS, although truncation of the outer disk by 
a close binary companion is also a viable possibility. The transition disks could be explained by a 
giant planet opening a gap in the disk, or a photoevaporative wind evacuating the inner disk. The 
latter mechanism requires a very low accretion rate and can be ruled out for the transition disk in 
our sample that is still accreting. 

None of the groups of evolved disks appear systematically older than the rest of the cluster, 
although the small numbers make any difference harder to evaluate. Initial conditions in the disk 
and surrounding envelope could influence how quickly disks evolve. The initial angular momentum 
or mass of the disk could also have an influence. 

This publication makes use of data products from the Two Micron All Sky Survey, which is 
a joint project of the University of Massachusetts and the Infrared Processing and Analysis Cen- 
ter/California Institute of Technology, funded by the National Aeronautics and Space Administra- 
tion and the National Science Foundation. This work is based in part on observations made with 
IRAC, aboard the Spitzer Space Telescope, and is supported by contract #1255094 for SSC/Caltech 
to the University if Arizona. The Second Palomar Observatory Sky Survey (POSS-II) was made 
by the California Institute of Technology with funds from the National Science Foundation, the 
National Geographic Society, the Sloan Foundation, the Samuel Oschin Foundation, and the East- 
man Kodak Corporation. We would also like to thank Tom Megeath for providing the IRAC data 
and George Rieke for comments that helped improve the paper. 
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Fig. 1 — POSS blue plate centered at a 2 ooo = 05 / M6 m 47 s <5 2 ooo = +00°09'50" spanning one degree 
on each side overlayed with the positions of cluster members. 
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Fig. 2. — SDSS color-magnitude diagram for those stars with optical photometry (dots) and those 
stars we classi fy as likely memb ers (squares). Most of the likely members lie above the main 
sequence locus (jSiess et al.ll2000l ). as expected for these young stars. 
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Fig. 3. — NIR color-magnitude diagram for all sources observed toward NGC 2068/71 (do ts) and 
all members (squares). The dashed line is the 2 Myr isochrone from Siess et al. (j200ol ) at the 
distance of NGC 2068/71. The solid line is a reddening vector extending from a 2 Myr MO star at 
the distance of NGC 2068/71, with every 5 Ay marked with an plus sign. 
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Fig. 4. — NIR color-color diagram of th e clusters members. Dashed line shows the m ain sequence 



locus from iKenvon &: Hartmannl (|1995l ) and the dotted line is the CTTS locus from iMever et al 



(|1997l ). The solid line is the reddening vector extending from an M6 star with every Ay = 5 
marked with a plus sign. The one star with anamolously blue colors is object 458 which may be 
contaminated by scattered light in the near-infrared. 




Fig. 5. — HR diagram for the likely members of NGC 2068/71. Different symbols represent dif- 
ferent types of disks. Blue symbols are CTTS and red symbols are WTTS. See text for definition 
of strong/IRAC-weak/MIPS-weak/no disk as well as description of derivation of luminosity and 
effective temperature. Dotted lines show mass tracks for a 0.3, 0.5, 1, 2 Mq star. A typical error 
bar is shown in the upper right corner. 
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Fig. 6. — Radial velocity histogram for the likely members. Dashed line shows the median velocity 
of likely members at 27.1 km/sec. Hashed histogram shows the CTTS in the cluster, which make 
up most of the likely cluster members around this median velocity. 
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Fig. 7. — Rotational velocity as a fraction of the break-up velocity (vbreakup = (GM / R) 1 / 2 ) versus 
mass. Upper limits are marked with arrows. Red symbols indicate stars with disks based on the 
presence of infrared excess and blue symbols represent stars without disks. 
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Fig. 8. — Sample hectochelle spectra focusing on the Ha profiles. The stars include a strong disk 
(63), a MlPS-weak disk (581), an IRAC-weak disk (515), a transition disk (281), the five stars 
with repeat hectochelle spectra (326, 458, 177, 1116 and 618) and one star without a disk (141). 
Stars #177, 141 demonstrate the narrow, symmetric profiles due to chromospheric emission while 
the rest show broad, asymmetric features due to accretion. Object ID numbers are included in 
each plot, along with a vertical line marking the line center, and the spectra of stars with multiple 
observations have been shifted vertically for easier comparison. 
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Fig. 9. — Equivalent width of the Ha line versus the full width at 10% maximum height for the 
Ha line for those likely members with high-resolution spectra. Vertical dashed line is located at 
250 km/sec, our boundary for actively accreting stars. Different symbols represent different types 
of disks while filled symbols are stars with a measurable amount of veiling. Stars with multiple 
spectra have the two observations connected by a line 
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Fig. 10. — Veiling, accretion luminosity and accretion rate as a function of mass. Upper limits are 
marked with arrows. See text for description of conversion from veiling to accretion rate. 
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Fig. 11. — Slope of IRAC SED versus effective temperature for the likely members. CTTS are 
marked with a square, WTTS are marked with a triangle. Dashed lines mark boundaries between 
stars with strong or MlPS-weak disks, stars with IRAC-weak disks and stars without disks. 
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Fig. 12. — IRAC color-color diagram for those stars in our spectroscopic sample. Size of the symbol 
indicates strength of Ha, as measured by its equivalent width, while the color of the symbol scales 
with its T e ff. The stars with the smallest infrared excess tend to have the smallest Ha EW, while 
the stars with the largest infrared excess in both colors have the largest Ha EW. 
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Fig. 13. — IRAC+MIPS color-color diagram for those stars in our spectroscopic sample. Size of 
the symbol indicates strength of Ha, while color indicates T e tf, as in figure [T2l 
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Fig. 14. — SED slope for the IRAC bands versus the SED slope from 8-24//m. The size of the 
points has been scaled to the size of the Ha equivalent width, while color scales with T e tt as in 
figure fl2l The dashed lines delineates between strong disks, IRAC weak disks and stars with no 
disk. The MIPS weak sources are also evident around cumips ~ — 1> an d the transition disks have 
large olmips an d small oirac- The SEDs have not been dereddened. 
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Fig. 15. — SEDs of stars exhib iting different le yels of infrared exces s emission. SEDs have been 
dereddened using Mathis (jl99ol ) for A < 3/iin and Flaherty et al. ( 2007 ) for A > 3/um and normalized 
to the H band flux. All of these stars have spectral type M1-M2. 
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Fig. 16. — Extinction map (axes in degrees) with class I sources from Muzerolle et al. (in prep) 
and the stellar members identified in this paper. The circle in the lower left shows the size of the 
gaussian kernel used to smooth the extinction map. The class I sources are located near the peaks 
in the extinction while the members follow th e broader distribution of dust. Contours correspond 
to Ay = 1, 2, 3, 4, 5 using the extinction law of lRieke Lebofskvl (119851 ) 
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a The source is not resolved from a nearby star in the J band and only an upper limit is reported. 

I 

CO 

Note. — SDSS u band data is not reported here because of the red leak in the filter. 



Table 2. Cluster members identified from Low- Resolution Spectra 
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CTTS 


288 


05 


45 53.60 


00 22 42.10 


K7 


K7 


2.4 


0.28 


WTTS 


^99 


05 


45 56.20 


UU OO 1U.OO 


M<\ MR 


HC r. 
ivu . o 


o.O 


n a 


VV -L -L O 


326 


05 


45 56.31 


00 07 08.58 


K7-M0 


K9 


6.2 


0.68 


CTTS 


344 


05 


45 57.38 


00 20 22.20 


K6-M0 


K7 


35.3 


0.5 


CTTS 


351+ 


05 


45 57.93 


00 02 48.55 


M3-M5 


M4 


75 


0.4 


CTTS 


373 


05 


46 00.18 


00 03 07.06 


M4 


M4 


27.2 


0.36 


CTTS 


416 


05 


46 04.58 


00 00 38.16 


M2 


M2 


7.4 


0.2 


WTTS 


418 a + 


05 


46 04.64 


00 04 58.15 


K7-M1 


K9 


126 


0.26 


CTTS 


458 


05 


46 07.89 


-00 11 56.87 


K2-K3 


K3 


3.2 


0.58 


CTTS 


480 


05 


46 09.27 


00 13 32.63 


M0-M1 


Ml 


10.8 


0.48 


CTTS 


515 


05 


46 11.86 


00 32 25.91 


K1-K2 


K2 


3.8 


0.7 


CTTS? 


543 


05 


46 14.48 


00 20 24.36 


M3-M4 


M4 


3.6 


<0.35 


WTTS 


571 


05 


46 18.30 


00 06 57.85 


K1-K3 


Kl 


34.7 


0.4 


CTTS 


581 


05 


46 18.89 


-00 05 38.11 


K3-K4 


K4 


3.9 


0.73 


CTTS 


584t 


05 


46 19.06 


00 03 29.59 


K4-K6 


K5 


44.5 


0.44 


CTTS 


590 


05 


46 19.47 


-00 05 20.00 


K2-K3 


K2.5 


28.5 


0.5 


CTTS 


601 


05 


46 20.88 


00 08 09.42 


K6-M1 


K8 


8.7 


0.7 


CTTS 



Table 2 — Continued 



Object ID 


RA (J2000) 


Dec (J2000) 


Spectral Type Range 


Adopted Spectral Type 


HaEW 

(A) 


Li EW 

(A) 


CTTS/WTTS 


618 


05 46 22.44 


-00 08 52.62 


K1-K2 


Kl 


32 


0.55 


CTTS 


628 


05 46 22.99 


00 04 26.44 


M1-M3 


M2 


104 


0.5 


CTTS 


656 


05 46 25.89 


00 09 31.97 


M1-M2 


M2 


32.4 


<0.4 


CTTS 


663 


05 46 26.65 


00 31 07.50 


K4-K8 


K6 


4.9 


0.6 


WTTS 


677 


05 46 27.83 


00 05 48.37 


M2-M3 


M2.5 


3 


0.3 


WTTS 


689 


05 46 29.00 


00 29 07.19 


M1-M2 


Ml. 5 


4.8 


0.5 


WTTS 


699 


05 46 30.06 


00 12 09.68 


M1-M2 


Ml. 5 


15.4 


0.5 


CTTS 


712t 


05 46 31.05 


00 25 33.37 


continuum 


- 


92.3 


0.25 


CTTS 


730 


05 46 33.28 


00 02 51.90 


K5-M0 


K7 


24.4 


0.6 


CTTS 


739 b 


05 46 34.54 


00 06 43.45 


K3-K5 


K4 


5.5 


0.5 


CTTS 


741 


05 46 34.90 


00 04 20.68 


M2-M3 


M2.5 


5 


0.2 


WTTS 


760 


05 46 37.06 


00 01 21.79 


K4-K7 


K5.5 


61 


0.4 


CTTS 


782 


05 46 38.57 


00 22 05.99 


M0-M1 


Ml 


5.8 


0.2 


WTTS 


802 


05 46 40.18 


00 05 01.86 


M0-M1 


Ml 


7.8 


0.8 


WTTS 


813* 


05 46 40.77 


00 27 22.50 


M2-M4 


M3 


166 


0.1 


CTTS 


843 


05 46 43.85 


00 15 32.36 


M2 


M2 


4.9 


0.6 


WTTS 


849 c 


05 46 44.09 


00 18 03.17 


K5-K7 


K6 


7.1 


0.7 


CTTS 


856 


05 46 44.84 


00 16 59.77 


M1-M2 


Ml 


14 


0.4 


CTTS 


878 f 


05 46 46.87 


00 09 07.63 


M0 


M0 


45 


0.4 


CTTS 


934 c 


05 46 51.48 


00 19 21.32 


M1-M2 


Ml. 5 


25.1 


0.3 


CTTS 


941 


05 46 52.41 


00 20 01.68 


K7-M1 


K9 


41 


0.5 


CTTS 


984 c 


05 46 56.54 


00 20 52.91 


K2-K4 


K3 


2.3 


0.7 


WTTS 


997+ 


05 46 58.03 


00 14 27.82 


M3-M4 


M3 


19 


0.45 


CTTS 


998 


05 46 58.13 


00 05 38.15 


K7-M0 


M0 


2 


0.73 


WTTS 


1062 


05 47 03.97 


00 11 14.35 


M1-M2 


M2 


42 


0.45 


CTTS 


1078 c 


05 47 04.94 


00 18 31.64 


M0-M1 


MO. 5 


2.3 


0.4 


WTTS 


1089 


05 47 05.34 


00 28 46.13 


M1-M2 


Ml. 5 


6 


0.3 


WTTS 


1099+ 


05 47 06.00 


00 32 08.48 


K0-K3 


KO 


20 


0.7 


CTTS 



Table 2 — Continued 



UDject ijj 


1\A [JZUDU ) 


uec [o zuuu ) 


Spectral Type Range 


Adopted Spectral Type 


(A) 


(A) 




1116 


05 47 06.96 


00 00 47.74 


K3-K6 


K4.5 


28.4 


0.5 


CTTS 


1117 


05 47 06.99 


00 31 55.96 


M0-M2 


Ml 


63 


0.3 


CTTS 


1119 c 


05 47 07.26 


00 19 32.23 


K7-M0 


MO 


48 


0.3 


CTTS 


1134 


05 47 08.69 


00 00 14.04 


M2-M4 


M3 


9.4 


0.42 


WTTS 


1171 


05 47 10.89 


00 32 05.96 


M1-M2 


M1.5 


8.4 


0.57 


WTTS 


1173 c 


05 47 10.98 


00 19 14.81 


G4-G8 


G6 


17 


0.6 


CTTS 


1211+ 


05 47 13.85 


00 00 17.06 


M0-M2 


Ml 


30.4 


0.7 


CTTS 


1212+ 


05 47 14.11 


00 09 07.34 


continuum 




78 


0.06 


CTTS 


1247 


05 47 16.58 


-00 00 56.38 


M4-M5 


M4 


10.3 


0.36 


WTTS 


1259+ 


05 47 17.16 


00 18 24.59 


M0-M2 


Ml 


65.5 


0.3 


CTTS 


1262 


05 47 17.28 


00 38 21.37 


K5-M0 


K7 


8.3 


0.6 


CTTS 


1284+ 


05 47 19.72 


00 01 21.83 


M5-M6 


M5 


300 


0.47 


CTTS 


1333 


05 47 25.05 


00 31 04.94 


K6-K7 


K7 


1.3 


0.24 


WTTS 


1643 


05 47 57.53 


00 11 31.24 


M3-M5 


M4 


9.5 


0.4 


WTTS 



f [O I] 6300 A emission detected in optical spectra 



a Variable in the optical as observed bv lBricefio et al.l (|2005l ). They measured a spectral type of K6, and Ha and Li EWs of 98 
and 0.4A respectively. 



'Binary from lPadgett. Strom fc Ghed (|l997h . 



c X-ray emission detected by ISkinner et al.l (j2007l ) 
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Table 3. Stellar Parameters for Likely Members 



Object ID 


Luminosity 


T eff 


Mass 


Radius 


A v 


Accretion Rate 


CHIRAC 


Disk? 




Lq 


K 


M Q 


Rq 




10- s M & yr- 1 






1 7 


n ooq 
U.zyo 


Q/1 70 

04 / U 


u.oo 


1 ^A 
1 .04 


1 

i.y 




O fi^ 
-Z.00 


no 


00 

zy 


U.UOo 


oUDU 


n it; 
U. 10 


u.yo 


o 
U.U 




o «a 

-z.04 


no 


Oo 


U.o4o 


Q/l 70 
o4 / U 


n Q/i 
U.o4 


1 RR 
1 .00 


1 1 
1.1 


O 70 

u. < y 


1 OQ 
-l.ZO 


strong 


1/11 
141 


U.014 


OoOU 


o c;7 
u.o / 


1 RA 
1 .04 


O Q 

z.o 


u 


O 70 
-z. 1 U 


no 


1 7/1 
1(4 


n 077 
U.z I I 


o4 / U 


U.o4 


1/18 
1 .4o 


1 

l.D 




77 
-z. / / 


no 


I 77 

I I i 


1 0^1 


a ^oo 


l.OU 


1 RK 
1 .00 


1 & 
l.o 


n 
u 




-Z.OO 


Xransition 


on/1 

zU4 


U.40O 


Q700 
o / ZU 


O /I 7 
U.4 / 


1 £7 
1 .0 ( 


1 1 
1.1 


O 

U.OO 


1 QO 
-l.oz 


strong 


000 

zzy 


n 077 
U.U ( < 


QQ70 


n o^i 
U.zO 


U.oo 


1 O 
l.z 




O 70 

-z. / y 


no 


281 


0.235 


3720 


0.47 


1.19 


i.i 


0.24 


-2.39 


Transition 


285 


0.348 


3580 


0.39 


1.56 


1.5 


0.87 


-1.12 


strong 


288 


0.230 


4060 


0.79 


0.99 


0.7 





-2.71 


no 


322 


0.214 


3240 


0.25 


1.49 


1.5 


- 


-2.61 


no 


326 


0.687 


3990 


0.68 


1.77 


1.6 


<0.47 


-1.79 


IRAC-weak 


344 


0.745 


3990 


0.69 


1.84 


1.6 


1.79 


-1.52 


strong 


351 


0.091 


3370 


0.28 


0.90 


0.8 


- 


-1.76 


IRAC-weak 


373 


0.496 


3370 


0.30 


2.10 


1.5 


1.65 


-1.15 


strong 


416 


0.319 


3580 


0.39 


1.49 


3.1 


- 


-1.45 


strong 


418 


0.895 


3990 


0.67 


2.01 


2.2 


2.27 


-1.03 


strong 


458 a 


4.679 


4730 


2.07 


3.33 


0.0 b 


<7.11 c 


-1.91 


IRAC-weak 


480 


0.807 


3720 


0.47 


2.20 


1.4 


1.77 


-1.13 


strong 


010 


Z.ZoU 


a onn 
4yuu 


1.00 


9 1 Q 

Z.lo 


i q 
l.o 


^-1 OO 

<i.uy 


9 9Q 

-z.zo 




04o 


U.11D 


oo / u 


U.oU 


1 .Ul 


l.o 




1 Q9 

-i.yz 


HVrVV^y-WedK 


0(1 


9 /1H7 
Z.4U 1 


4yUU 


1 RK 
l.DO 


9 10 

z.iy 


A 7 
4. / 


^-i 1 n 

<^l. 1U 


-U. I 


strong 


Ool 


Q £7£ 
0.0(0 


4oyu 


1 KR 
l.OD 


O.U4 


Q Q 
O.O 


07 
<z.o / 


1 9p; 

-l.ZO 


niTpq ...„„i, 
ivii.ro- weaK 


584 


3.280 


4350 


1.12 


3.244 


4.5 


7.83 


-1.06 


strong 


590 


3.525 


4815 


1.85 


2.74 


2.6 


3.22 


-0.72 


strong 


601 


0.963 


3955 


0.64 


2.13 


4.3 




-0.75 


strong 


618 


2.053 


5080 


1.26 


1.88 


2.6 


<1.06 


-1.67 


IRAC-weak 


628 


0.391 


3580 


0.39 


1.65 


3.4 




-1.02 


strong 


656 


0.352 


3580 


0.39 


1.57 


3.1 




-1.24 


strong 


663 


1.419 


4205 


0.89 


2.28 


4.4 




-1.27 


MlPS-weak 


677 


0.643 


3525 


0.37 


2.19 


4.1 




-2.64 


no 


689 


0.721 


3650 


0.43 


2.16 


2.2 





-2.66 


no 


699 


0.899 


3650 


0.43 


2.41 


4.6 




-1.21 


strong 


730 


1.280 


4025 


0.72 


2.37 


3.9 




-1.06 


strong 


739 


3.182 


4590 


1.63 


2.87 


6.2 




-1.16 


MlPS-weak 
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Table 3 — Continued 



Ohiect TD 


T iiimirinQitv 


rp 

1 eff 


Mass 


T? e A (\ \ 11 Q 






a IRAC 


Disk? 






K 




.n,0 




1U JW0J/I 






741 


0.625 


3525 


0.37 


2.16 


4.0 


_ 


-1.82 


IRAC-weak 


760 


1.655 


4277 


1.01 


2.38 


3.5 


5.13 


-1.13 


strong 


782 


0.491 


3720 


0.47 


1.72 


4.0 


_ 


-1.39 


MlPS-weak 


802 


0.552 


3720 


0.47 


1.82 


2.7 




-2.19 


IRAC-weak? 


813 


0.335 


3470 


0.34 


1.63 


3.2 


- 


-1.28 


strong 


843 


0.399 


3580 


0.39 


1.67 


2.4 


- 


-1.56 


strong 


849 


2.006 


4205 


0.90 


2.71 


2.9 


<1.79 


-1.26 


MlPS-weak 


856 


0.618 


3720 


0.47 


1.93 


4.2 


- 


-1.83 


IRAC-weak 


878 


0.605 


3850 


0.57 


1.78 


3.2 


1.09 


-0.73 


strong 


934 


0.422 


3650 


0.42 


1.65 


3.2 


- 


-1.09 


strong 


941 


0.505 


3920 


0.62 


1.57 


2.8 


1.54 


-0.96 


strong 


984 


2.094 


4730 


1.89 


2.19 


4.5 


<1.78 


-2.44 


IRAC-weak? 


997 


0.315 


3470 


0.34 


1.58 


3.4 




-0.93 


strong 


998 


3.269 


o o c n 
3850 


0.56 


A 1 O 

4.13 


3. I 


<5.44 


-1.79 


IKAU-weak 


1062 


0.441 


3580 


0.39 


1.76 


2.0 


1.94 


-0.91 


strong 


1078 


1.973 


3785 


0.51 


3.32 


4.2 


_ 


-2.24 


IRAC-weak? 


1089 


0.299 


3650 


0.42 


1.39 


3.4 


_ 


-2.64 


no 


1099 


2.361 


4900 


1.64 


2.17 


5.5 


_ 


-1.11 


strong 


1116 


1.532 


4470 


1.29 


2.10 


1.8 


<0.81 


-0.97 


strong 


1117 


0.362 


3720 


0.47 


1.47 


3.8 




-0.56 


strong 


1119 


7.294 


3850 


0.56 


6.18 


3.8 


<18.2 


-0.80 


MlPS-weak 


1134 


0.667 


3470 


0.35 


2.30 


1.2 





-2.62 


no 


1171 


0.519 


3720 


0.47 


1.74 


2.4 




-1.28 


MlPS-weak 


1173 


9.082 


5700 


1.84 


3.14 


5.1 


no stand 


-1.13 


strong 


1211 


1.302 


3720 


0.47 


2.80 


1.6 


7.75 


-1.51 


strong 


1247 


0.253 


3370 


0.30 


1.50 


2.5 




-1.29 


MlPS-weak 


1259 


0.336 


3720 


0.47 


1.42 


3.3 




-1.37 


strong 


1262 


2.798 


4025 


0.70 


3.50 


4.5 


8.60 


-1.31 


strong 


1284 


0.108 


3240 


0.22 


1.06 


1.4 




-1.10 


strong 


1333 


0.216 


4060 


0.79 


0.96 


0.5 





-2.71 


no 


1643 


0.348 


3370 


0.30 


1.76 


1.9 




-2.62 


no 



a Photometry may be contaminated by background source with very cool SED, possible a Class 0/1 
source. 

b With no available optical photometry, and extended emission contaminating the 2MASS photom- 



etry, an accurate extinction could not be derived and we assume a lower limit of 



C R band magnitude from ISiess et al.l (J2000) 2 Myr isochrone used instead of observed R band 
magnitude to derive accretion rate from veiling. 



Note. — A dash in the accretion rate column indicates no high-resolution spectra was available, a 
is used for those stars with high-resolution spectra, but the Ha profile indicated there was no active 
accretion. Disk descriptors: strong {aiRAC> —1-6, [8.0]-[24]> 2.4), MlPS-weak (amAC< —1-6, 
[8.0]-[24]< 2.4), IRAC-weak( — 1.6 >ain.AC> — 2.5), gap (small short wanvelength excess, large [24] 
excess), no disk (a/KAC< —2.5). 
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Table 4. High-Resolution Spectra of Likely Members 



Object ID 


Radial Velocity 
(km/sec) 


Vsini 
(km/sec) 


T&D R a 


Ha FW 10% 

(km/sec) 


HaEW 
A 


Veiling 


63 


31.7±3.0 


<8 


3.7 


280 


19.1 


0.88 


141 


27.0±1.3 


<8 


12.0 


142 


1.7 


- 


177 


27.7±0.8 


<8 


23.6 


124(101) 


0.9(1.2) 


- 


204 


25.5±1.2 


<8 


11.5 


250 


25.4 


0.51 


281 


27.3±1.3 


<8 


16.7 


400 


32.4 


0.36 


285 


26.2±3.6 


14.8 


3.9 


600 


120.7 


0.94 


288 


30.2±1.3 


<8 


13.4 


110 


1.4 


- 


326 b 


27.0±1.1 


<8 


15.5 


170(170) 


5.9(5.3) 


- 


344 


21.6±1.3 


<8 


12.1 


495 


27.4 


0.76 


373 








290 


24.8 


2.2 


418 


27.0±2.1 


13.4 


7.4 


410 


60.5 


0.76 


458 


24.5±1.1 


29.7 


18.4 


280(320) 


2.7(3.4) 


<0.15 


480 


26.8±1.3 


<8 


12.5 


510 


7.9 


0.45 


515 b 


14.1±1.0 


24.3 


20.1 


410 


3.4 


<0.15 


571 


28.8±1.5 


15.7 


11.3 


483 


29.6 


<0.15 


581 


25.9±0.9 


18.3 


21.2 


480 


4.8 


<0.15 


584 








640 


42.7 


0.39 


590 


28.2±0.7 


<8 


23.1 


470 


14.9 


0.31 


618 


27.0±0.7 


26.8 


19.2 


530(520) 


46.8(28.3) 


<0.15 


689 


27.2±1.2 


<8 


10.9 


110 


3.7 




760 


28.5±1.8 


<8 


8.2 


535 


49.4 


1.28 


849 


27.3±1.6 


<8 


9.8 


340 


18.2 


<0.15 


878 


27.1±1.3 


<8 


12.6 


440 


25.5 


0.49 


941 


26.8±2.3 


<8 


6.5 


415 


52.4 


1.82 


984 


25.4±2.4 


49.2 


10.6 


280 


1.1 


<0.15 


998 


11.6±2.3 


40.7 


9.7 


280 


1.9 


<0.15 


1062 


26.9±2.9 


<8 


4.3 


430 


46.8 


3.18 


1116 


25.8±0.8 


<8 


21.9 


400(405) 


23.4(22.1) 


<0.15 


1119 


31.3±3.8 


55.8 


2.9 


450 


27.1 


<0.15 


1134 


27.1±1.5 


<8 


8.7 


160 


8.0 




1173 


32.7±2.3 


49.1 


11.4 


350 


9.7 
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Table 4 — Continued 



Object ID 


Radial Velocity 
(km/sec) 


Vsini 
(km/sec) 


T&D R a 


Ha FW 10% 

(km/sec) 


HaEW 

A 


Veiling 


1211 


25.9±2.5 


<8 


5.2 


490 


26.3 


1.97 


1262 


26.2±1.2 


<8 


14.1 


340 


5.4 


0.41 


1333 


28.3±1.1 


<8 


15.6 


90 


1.3 





Note. — FW at 10% values in parenthesis are from the second set of high-resolution 
spectra. Only objects 326,177,458,618 and 1116 have two sets of high-resolution spectra. 



a Tonry & Davis R value (jTonrv David Il979l ). The 90% confidence interval for v 
sini is approximately vsini/ (1+R). 

b Ha profile shows an inverse P Cygni shape indicating active accretion, despite the 
small FW 10%. 



Table 5. Fraction of Evolved Disks 



Cluster Age Strong disks MlPS-weak IRAC-weak 



K0-M1 



Taurus 


1 


94 ±24% (16/17) 




6 ±6% (1/17) 


NGC2068/71 


2 


66 ± 14% (23/35) 


16 ±7% (5/32) 


20 ± 8% (7/35) 


IC 348 


2-3 


80 ±20% (16/20) 


0% (0/17) 


20 ± 10% (4/20) 


cr Ori 


3 


75 ± 19% (15/20) 


19 ±11% (3/16) 


10 ± 7% (2/20) 


Tr 37 


4 


68 ± 9% (59/87) 


13 ± 5% (7/54) 


24 ±5% (21/87) 


NGC 2362 


5 


30 ± 17% (3/10) 




70 ± 26% (7/10) 



Note. — Statistics for MlPS-weak disks only include those stars detected at 
24/im. IRAC-weak disks denned as —2.56 <otiRAC< —1.8. 
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